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ABSTRACT 

Galactic neutral hydrogen (HI) within a few hundred parsecs of the Sun con- 
tains structure with an angular distribution that is similar to small-scale structure 
observed by the Wilkinson Microwave Anisotropy Probe ( WMAP). A total of 108 
associated pairs of associated HI and WMAP features have now been cataloged 
using HI data mapped in 2 km/s intervals and these pairs show a typical offset of 
0.°8. A large-scale statistical test for a direct association is carried out that casts 
little additional light on whether the these small offsets are merely coincidental 
or carry information. To pursue the issue further, the nature of several of the 
features within the foreground HI most closely associated with WMAP struc- 
ture are examined in detail and it is shown that the cross-correlation coefficient 
for well-matched pairs of structures is of order unity. It is shown that free-free 
emission from electrons in unresolved density enhancements in interstellar space 
could theoretically produce high-frequency radio continuum radiation at the lev- 
els observed by WMAP and that such emission will appear nearly fiat across the 
WMAP frequency range. Evidence for such structure in the interstellar medium 
already exists in the literature. Until higher angular resolution observations of 
the high-frequency continuum emission structure as well as the apparently associ- 
ated HI structure become available, it may be difficult to rule out the possibility 
that some if not all the small-scale structure usually attributed to the cosmic 
microwave background may have a galactic origin. 

Subject headings: Interstellar matter, neutral hydrogen, cosmology, WMAP 



1. Introduction 



The interpretation of the distribution of the small-scale structure observed by the 
Wilkinson Microwave Anisotropy Probe (WMAP), as epitomized by the summary prepared 
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from the five-frequency data in what is called the Internal Linear Combination (ILC) map, 
forms a cornerstone of modern cosmology. The ILC map has been presented by Hinshaw et 
al. (2007) and has gone through several iterations, referred to there. The key aspect of the 
ILC map is that the observed structure appears to be consistent with the existence of sound 
waves in the early moments of the universe, as summarized in the shape of the so-called 
acoustic spectrum. 

It was disturbing to the present author to discover that some of the small-scale struc- 
ture in the ILC data appeared to be closely associated with small-scale structure in the 
distribution of interstellar neutral hydrogen (HI) emission in the Galaxy, Vcrschuur (2007a, 
hereafter Paper 1, & 2007b). If these associations were to be regarded as anything other 
than mere coincidence they would imply that a previously unrecognized process occurring 
in interstellar space is capable of generating the structure in the high-frequency continuum 
emission (HFCE) observed by WMAP. That, in turn, would imply that the cosmological 
interpretation of the ILC structures would be affected. 

An immediate criticism of Paper 1 was produced by Land & Slosar (2007) who searched 
for one-to-one correspondence between HI and ILC features. However, that did httle to either 
confirm or invalidate our earlier discussions, which were based on the discovery that peaks 
in the two forms of emission are closely associated but generally offset by a small angle. This 
claim will be re-examined through a statistical test for a point-to-point correlation between 
the two forms of emission in narrow bands of HI emission. When HI structure is mapped with 
a velocity resolution of 3 km/s at intervals of 2 km/s, this produces area maps at 200 different 
velocities between —200 to -1-200 km/s. (All velocities are given with respect to the local 
standard of rest.) It is the apparent association between ILC structure and HI features in 
these velocity maps that reveals the possibility that there exist physical processes occurring 
in interstellar space capable of creating the observed HFCE. The question that must then be 
considered is whether there is any possibility that such processes can occur in the interstellar 
medium. It will be shown that unresolved structures containing excess electron density could, 
in principle, give rise to the HFCE at the levels observed by WMAP. Evidence for structures 
of the required angular scale, particle density and temperature already exists in the literature. 

The present analysis elaborates on the work reported in Paper 1 where it was pointed 
out that toward HI cloud MI there is a close association between ILC and HI structure 
as well as excess soft X-ray emission reported by Herbstmeier et al. (1995). MI is one of 
two HI features at anomalous velocities discovered by Mathewson (1967) that came to be 
named MI and Mil, the former at velocities around —110 km/s and the latter at about —80 
km/s. In the case of both MI and Mil, associated Ha emission has been reported by Tufte, 
Reynolds & Haffner (1998). In Paper 1 it was predicted that a similar relationship between 
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/LC, HI and X-ray structure found for MI would also be manifested in MIL The prediction 
is confirmed. 

In §3 the new statistical test is described. In §4 tables of HI-/LC associations in an 
area of sky below Mil as well in directions of the ten brightest ILC peaks in an area of sky 
encompassing all galactic longitudes between | h\ — 30° & 70° are presented. In §5 the average 
angular separation between HI and associated /LC peaks found in Paper 1 and in the body of 
this paper are discussed. In §6 the prediction made in Paper 1, that a relationship would be 
found between ILC structure, HI and soft X-ray emission for Mil is discussed. The results of 
Gaussian analysis of the profiles toward HI features MI & Mil are outlined in §7 in order to 
determine if something can be learned about the relationship of the narrow-band HI features 
to the ILC structure. In §8 a number of the very close associations between HI and ILC 
structure in the vicinity of both MI and Mil are considered and in §9 the results of cross- 
correlation calculations for several pairs of structures are presented. In §10 it is shown that it 
is theoretically possible that unresolved structure in the distribution of interstellar electrons 
can give rise the high-frequency emission at the same level as observed by WMAP. In §11 a 
variety of data pertaining to the presence of small-scale interstellar structure is considered 
in the light of the present work. Discussion is offered in §12 followed by conclusions in §13. 



2. Data 

The HI data used in this study were taken from the side-lobe-corrected, Leiden- Argentina- 
Bonn {LAB), All-Sky HI survey carried out with a beamwidth of 0.6° and bandwidths 1.3 
km/s (Kalberla et al. 2005). These data were used to produce contour maps in galactic lon- 
gitude and latitude (called /, h maps) of the total HI column density as well as the brightness 
temperature integrated over a 3 km/s velocity range plotted every 2 km/s in velocity over 
the relevant areas of interest. During this study it became apparent that while a 10 km/s 
interval was used in Paper 1, a great deal more information emerges when the emission in 
a 3 km/s wide band plotted at 2 km/s intervals is used. In addition, position-velocity plots 

or h,v maps) were produced when needed and in several cases individual HI emission 
profiles were extracted from the LAB database for the purpose of Gaussian analysis to be 
described. 

In this report we will use the summary WMAP data produced by Hinshaw et al. (2007) 
who combined the five-frequency data (in the range 23 to 94 GHz), after correcting for a 
number of effects including contributions from the galactic foreground, to produce the so- 
called Internal Linear Combination {ILC) map whose angular resolution was smoothed to 1°. 
Structure seen in this map is generally taken to represent conditions in the early universe. 
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In order to understand the causes of the apparent associations between certain ILC and 
galactic HI features, observational studies should focus on those directions where the HI 
structure is relatively uncomplicated. This rules out all directions close to the galactic disk. 
Therefore the nature of the HI- ILC relationship in the galactic latitude band bordered by 
30°and 70° in both galactic hemispheres has been studied, the lower limit chosen to avoid the 
galactic disk, the upper limit set for practical purposes to avoid getting too close to the poles 
given that we use the rectangular galactic coordinate system of the LAB survey. The bulk 
of the work reported here focusses on the northern galactic hemisphere data, particularly 
between I — 160° & 210°. Examination of the data shown by Hinshaw et al. (2007) shows 
that this area of sky was subjected to the smallest corrections for foreground effects. 

The first year ILC data made available to the author for the work initiated in Paper 1 
were used for the present analysis and were supplemented by data for all longitude between 
latitudes 30° & 70° for both galactic hemispheres. Note that at the high galactic latitudes to 
which this work is confined there are no significant differences in the morphology of the ILC 
structure between the 1st and 3rd year WMAP results. After all, very little if any correction 
for foreground continuum radiation (Hinshaw et al. 2007) was applied to the high galactic 
latitude data other than toward the spurs of galactic continuum emission, which lie outside 
the bounds of the regions considered in detail in this report. 



3. Large-scale statistical tests 

Given that nearly every (positive amplitude) ILC peak examined closely in Paper 1 and 
in what follows appears to be associated with an HI peak at some velocity, statistical tests to 
confirm or negate their significance may be moot. The point is that no amount of statistical 
testing can prove a negative. However, in order to determine if there is any evidence for a 
direct point-to-point correlation between the ILC and HI structure data, a simple statistical 
test was performed. On the one hand the ILC data are available as intensities as a function 
of ^ & b. Hence we have Ijlc {ijb). On the other hand the HI data are available as a function 
of velocity as well. Thus the HI data are available in a data cube containing the brightness 
temperature, T^ as a function of l,b, and v. Hence Tb (/,&,?;). The statistical test involves 
examining the product: 

P{v)^TB{l,b,v).IjLc{l,b) (1) 

at each of the velocities at which galactic emission HI over the area of interest is found. 
This product was calculated using HI data integrated over a 3 km/s band (2 channels) at 
intervals of 2 km/s from —200 to +50 km/s over the northern sky between the latitude 
limits of 30° & 70° and between —100 and +50 km/s for the southern sky data between 
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the same latitude limits. These limits were set by the extent of the HI emission in the two 
hemispheres. Furthermore, the data were divided into thirds for areas in the northern (N) 
and southern (S) hemispheres as follows: 

Target Area I = 60°to 180° (TN or TS) 

Comparison Area 1 / = 300°to 60° (CIN or CIS) 

Comparison Area 2 / = 180°to 300° (C2N or C2S) 

The appellation " Target Area" grew out of the work reported in Paper 1 which focused 
on that longitude range in the north. 

In addition to calculating the product, P{v), for a given area, comparison calculations 
were performed by taking the HI data for that area, say TN, and then overlaying it on the 
/LCdata for each of the two comparison areas, CIN and C2N, and repeating the calculations. 
Thus three sets of products are calculated for each of the areas hsted above: 

P{v)tn.tn = Tb{1, b, v)tn-Iilc{1, b)TN (2) 

-P(^^)rjv.cijv = Tb{1, b, v)TN-hLc{l, b)ciN (3) 
P{v)tn.C2N = Tb{1, b, v)tn-Iilc{1, b)c2N (4) 

Similarly, three sets of products were calculated for each of the Comparison Areas 1 & 
2. The entire process was repeated for the data in the southern galactic hemisphere. It is 
recognized that the HI and ILC data have different resolutions but that does not negate this 
search for apparent associations since the higher resolutions HI structures would appear to 
be encompassed by the lower resolution ILC data if the structures happened to overlap. 

If the distribution of the WMAP structure and the HI peaks are random with respect 
to one another, a histogram of the product, 

Piv)^TB{l,b,v).IiLcil,b), (5) 

should show a normal distribution. If, however, there is a significant correlation between 
the two forms of emission within some velocity range the histogram should show a relative 
excess above a normal distribution at those velocities. 

Each group of three histograms covering 120° in longitude and calculated every 2 km/s 
in velocity produced roughly 3 x 10^ data points. For display purposes the results were 
plotted in contour map form for the number (N) of occurrences of the product with a given 
value (P) as a function of velocity. Hence we display N(P,?;). For a given one third of the area 
studied (according to equations 2, 3 & 4), three contour maps were plotted. This produced 
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a total of 9 such maps for the northern hemisphere and another 9 for the southern sky data. 
In all, the full set of calculations involved 2.7 x 10^ individual data points. 

Two sets of values for the product were calculated for each data set. One considered 
the value of P in steps of 0.05 K.mK and another over a limited range centered on zero in 
steps of 0.01 K.mK in case those revealed more detailed structure. 

In Paper 1 it was claimed that there exists a relationship between ILC structure with 
positive values and small-scale HI structure. ILC structures that have positive values are 
deemed to represent directions of a small excess of high-frequency continuum radiation above 
the 2.738K cosmic microwave background and negative values represent slightly cooler areas, 
with the range from —0.22 to -1-0.22 mK. 

We have previously noted (and will discuss this further below) that there is a typical 
offset between HI and WMAP ILC peaks of order 0.°8. If this offset were larger than the 
typical width of the ILC and HI structures, no correlations should make their presence known 
in the map of the product, P(v). If, instead, there is overlap between the contours defining 
structure in the two forms of emission the product, P(v), should reveal the presence of such 
associations in certain velocity ranges. 

Figure 1 shows three histogram contour maps for the HI data for the Target Area in 
the north (TN) when correlated with the ILC data for Comparison Area 2 (at the left), 
with the ILC data for the Target Area (center), and when compared with the ILC data 
for Comparison Area 1 (right hand plot). The product calculations were binned in steps of 
0.01 K.mK in these plots. Asymmetries toward negative product values implies a correlation 
between the HI structure and negative ILC structure. These diagrams appear to reflect 
an excess of negative value signals in the ILC data in the northern strip of sky considered 
here. Significant cases of direct positional associations between the HI and ILC data in 
the center plot should show up as marked asymmetries in the value of the product over a 
range of velocities pertaining to the relevant HI. However, very little asymmetry is found in 
these plots bearing in mind that such small asymmetries as may be visible involve relatively 
few data points compared to the peaks in the histograms. The individual histograms are 
themselves seldom simple Gaussians in shape. This, too, is not surprising, since neither 
the distribution of large-scale distribution of HI and ILC structures over tens of degrees 
are not randomly distributed on the sky. This is evident from a visual examination of the 
available data. It is the small-scale associations that are of interest. For example, the shght 
asymmetries around —50 km/s in the left-hand and right-hand plots indicate that in some 
areas of sky areas of excess HI structure happens to overlap areas of predominantly cither 
positive or negative patches of ILC structure. Thus taken together the data do not indicate 
significant overlap in position of the HI and ILC features at any favored velocity. 
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Land & Slosar (2007) have performed a statistical test using HI data integrated in 
velocity bands 10 km/s wide over the whole sky and performed a search for direct correlations 
between ILC and HI structure in these velocity intervals. However, any associations that 
might exist between the two forms of emission would not be expected to occur at a specific 
velocity over large areas of sky. 

The three plots in Fig. 1 and the other 15 similar plots that were produced for the 
northern and southern galactic hemispheres between |6| = 30° & 70° show that convincing 
evidence for one-to-one associations is absent. This docs not contradict the presence of 
associations that manifest as small angular offsets between HI and ILC peaks. Those are 
the ones that should be the focus of any further study. The point is that close associations 
must first be identified by actually examining the relationship between the ILC and HI area 
maps and then calculating a cross-correlation coefficient on a case-by-case basis, taking into 
account that the associated peaks are usually offset by a small amount and that the axis 
joining the offset structures do not share a common position angle over the sky. 

4. Newly Identified HI-/LC Pairs 

In this section a number of HI-/LC pairs will be cataloged in an area where the apparent 
association between the two types of structure is very obvious, even to the naked eye when 
examining the relevant HI and ILC data, for example the all-sky HI map found in the LAB 
web site (see astro.uni-bonn.de) and the widely publicized ILC map (see map.gsfc.nasa.gov). 

4.1. A Table of Associations in the Mil Area 

HI feature Mil is located at {l,h) = (185°, 63°) and its relationship to ILC structure was 
studied in the context of a larger area bounded by / = 180° & 210°, & = 50° & 70°. Figure 
2a shows the total HI content for this area with the ILC contours overlain. (In this and 
subsequent figures the HI data will be shown in color in the on-line version and as inverted 
gray-scale in the print edition.) Figure 2b shows the same HI data in contour map form 
for visual comparison with Fig. 2a. Obviously one cannot overlay two contour maps of 
this complexity in one figure and expect to learn anything by visual inspection alone. The 
learning comes from taking the two contour maps and examining similarities and differences 
between them in a "blink comparison" mode. Some readers may prefer that such comparison 
be undertaken by computer, but one can only program the computer to find what one is 
looking for, and one cannot determine what one is looking for without first looking at the 
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data! 



4.2. A Trio of Associated Structures neeir Mil 

In Fig. 2 only a few associations between HI and ILC peaks are obvious. This impression 
changes dramatically when the HI data in 3 km/s wide bands are plotted at 2 km/s intervals. 
For example, Fig. 3a shows a very bright (relative to its environment) HI feature at {l,b) — 
(201. °5, 56. °5) in the velocity range —50 to —48 km/s which is all but invisible in the map 
of total column density (Fig. 2). This HI peak is associated with an ILC peak observed 
by WMAP at {l,b) = (201. °2, 56.°6) which is listed as #79 in Table 1, below. Figure 3b 
shows an elongated peak at {l,b) = (203°, 55°) found in the HI integrated between —40 & 
—38 km/s. Fig. 3c shows the HI brightness integrated between +1 and +3 km/s and it is 
associated with the prominent JLC peak at {l,b) = (206. °4, 55. °0) listed as #74 in Table 1. 
In Fig. 3d the HI integrated brightness distribution in these three velocity ranges are plotted 
together and shows how closely the HI and ILC structures are correlated even if the HI is 
found at distinctly different velocities. Note that the shght offsets in longitude between the 
small angular diameter peaks in the HI and ILC data seen in Fig. 3 are 0.°3, or half the 
HI beam width. Further details in the HI distribution in narrow velocity intervals will be 
discussed in §6, 7 & 8 below. 

4.3. A Table of Associations 

The 34 brightest peaks in the ILC contour map have amplitudes >0.063 mK and their 
positions and amplitudes are given in Table 1 together with the properties of the associated 
HI peaks found in the narrow-band maps. Column 1 assigns an identification number for 
the ILC peaks in the Mil region while continuing the numbering scheme begun in Paper 1, 
Table 1. The longitude and latitude of each ILC peak is indicated in columns 2 & 3 and the 
amplitude of the peak in mK in Column 4. The center velocity of the 3 km/s wide band in 
which the associated HI feature was recognized is given in Column 5. Columns 6 & 7 give 
the longitude and latitude of the relevant HI feature and its peak amplitude in K.km/s in 
Column 8. The angular offset in arcdegrees between associated pairs of peaks is given in 
Column 9. 

In numerical terms 31 of the 34 ILC peaks in Table 1 reveal a closely spaced HI peak 
at some velocity when using these HI data plotted in 2 km/s intervals. (This contrasts with 
only 7 of the 34 that appear to be associated in the HI map of total column density. Fig. 2.) 
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In all, 48 distinctly different HI features in the narrow-band HI maps appear to be associated 
with these 31 ILCpeaks. This is consistent with what was reported in Paper 1, which shows 
that in many cases HI at more than one velocity is apparently involved in the production of 
the HFCE, possibly at the interface between interacting HI features. The average separation 
between paired set of Hl-ILC peaks listed in Table 1 is 0.°67 ± 0.°35, which happens to be 
closely equal the beamwidth of the LAB HI survey. 

The individual HI maps for this area between them contain 71 HI peaks defined as 
sets of closed contours that can be followed over adjacent velocity intervals, many of them 
showing small velocity gradients with changing position. About 23 of these peaks, or 32%, 
identified in the individual maps do not appear to have associated positive amplitude ILC 
peaks. This compares with 6% of the ILC peaks that do not show associated HI features. 

A comparison was also carried out for the relationship between HI peaks and negative 
ILC peaks with very similar results. The area of Fig. 2 contains 30 distinct minima with 
amphtudes < 0.056 mK and 18 of those appears to be associated with HI peaks. The average 
separation is 0.°90 ± 0.°34. Given the high density of structures in the ILC data for this 
region, with clear minima separating the peaks, it is to be expected that in general the 
minima would tend to mimic what is found for the maxima as regards their offsets with 
respect to HI peaks. One way to become certain that the associations between the positive 
value ILC peaks and HI that are significant is to determine if they reveal anything new 
about interstellar physics. Better still would be to obtain higher resolution observation of 
the HFCE and the HI structure to examine their apparent relationship more closely. 



4.4. The Ten Brightest ILC peaks 

The ten brightest ILC features in the area of sky for all galactic longitudes between b 
= 30° & 70° in both the northern and southern galactic hemisphere were also examined to 
determine if they were associated with HI structure. All 10 sources, nine of which are found 
in the southern sky, arc listed in Table 2 where the columns refer to the same properties as 
in Table 1. The numbering in Column 1 continues from the numbering used in Table 1. The 
velocities of the associated HI were again derived from examination of maps made at 2 km/s 
intervals. Nine of the ten brightest ILC peaks can readily be related to the presence of HI 
and most of those involve HI peaks at more than one velocity, which is similar to what is 
found in Table 1. 
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5. Average Sepeirations of Associated Pairs 

The preliminary studies of possible associations between HI & ILC structures discussed 
here and in Paper 1 so far includes 108 close pairs. (An overview of the southern sky data 
for which we have prepared maps to continue this work reveal at least another 100 obvious 
close associations. Details are deferred to a later paper.) Table 3 summarizes the average 
angular separation between the ILC and HI peaks in each pair taken from three independent 
data sets considered so far. The first entry is obtained from Table 1 in Paper 1. The second 
entry is from Table 1 above and the third entry is from Table 2 above. 

The striking fact that emerges from the data in Table 3 is that the average apparent 
angular separations between HI peaks associated with WMAP ILC peaks are all of order 
0.°8. This may be related to the method in which the associations were identified, which 
required that the HI contours defining a peak exhibit a morphology similar to the ILC 

feature for contours >0.02 mK. For cases where the separations are larger than about 2°, 
confusion quickly prevents an association from being recognized. At that separation the ILC 
amplitudes tend to become negative. 

At this point we reiterate concerns that are raised by the overall situation summarized 
in Table 3. In total, of the 108 ILC peaks 102, or 94%. show associated HI peaks. The 
converse is not true in that not all HI peaks show assocciated ILC structure. These cataloged 
associations involve HI at 152 distinct center velocities, which suggests that it could be the 
interaction between HI at different velocities that gives rise to the weak HFCE observed by 
WMAP. Given that the HI is pervasive and extremely patchy, this abundance of associations 
is not surprising and could obviously be fortuitous. The issue then becomes one of deciding 
whether it is possible that a previously unrecognized process in interstellar space could 
produce the HFCE observed by WMAP (see §10). 

6. A Close Look at the HI-/LC Association in the MI & Mil Areas 

In Paper 1, Fig. 6 the relationship between ILC, HI and excess soft-ray emission was 
pubhshed for MI. Those data are reproduced here in Fig. 4a in which the shaded pixels 
representing the X-ray data taken from Herbstmeier et al. (1995) are plotted on a map of 
total HI content between velocities of —140 & —80 km/s with the ILC contours overlain. 
The nature of the HI emission underlying the shaded pixels will be discussed in §7 below 
because here the brightest HI features are hidden by the X-ray pixel structure. In Paper 1 
it was predicted that something similar would be found for Mil but the data that allowed 
that assertion to be tested only came to hand after that paper was submitted. A similar plot 
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for Mil is shown in Fig. 4b where the HI data were integrated between —125 & —55 km/s. 
[Note that in these two areas very httle HI emission is seen at low velocities, as reported by 
Verschuur & Schmelz (2010a).] Again it appears that all three forms of emission are related 
for Mil and while in MI all three are closely associated in position for Mil the peaks for each 
type of radiation are slightly offset from the other two. If these close associations are due to 
anything other than pure chance then their relative location with respect to each other must 
surely provide information about the underlying physical processes involved in producing 
the HFCE and soft X-rays associated with the HI structure. 



7. Gaussian Analysis of HI Profiles Toward MI &; Mil 

In order to more closely explore the possible relationship between HI and ILC structures 
revealed in the data, HI emission profiles in several areas of interest were decomposed into 
Gaussian components and the column densities of families of component line widths sepa- 
rately mapped on the sky. This Gaussian decomposition was done by paying attention to the 
existence of underlying component 34 km/s wide that appears to be present in most if not all 
directions observed to date in the northern galactic hemisphere as reported in a preliminary 
study by Verschuur (2007b). This broad component appears to be separately present in each 
of the so-called high-, intermediate-, and low-velocity regimes, hereafter referred to as HV, 
IV and LV. This component is described by Verschuur & Pcratt (1999) and Verschuur (2004) 
and is readily identified in directions of the simplest HI profiles, while Verschuur (2007b) 
showed that in certain directions where the HI is least complex the pervasive nature of this 
component is most obvious. Therefore, in the present study a Gaussian component with this 
width was fit to each HI profile in each velocity regime while the Gaussian fitting algorithm 
(Verschuur, 2004) solved for the other components present. 

Note that it is possible that a Gaussian decomposition of a single profile may give 
ambiguous results. However, when Gaussian analysis over an area is carried out, at least for 
relatively simple profiles, a remarkable coherence emerges when the results are considered 
as a whole. Verschuur (2004) and references therein focused on relatively simple profiles 
to avoid ambiguity. In this respect the profiles shown in Fig. 4 are in fact simple in that 
the three main peaks are well separated in velocity. Thus each of those peaks is in itself 
simple, requiring at most three components to obtain a fit. Severe problems related to 
ambiguity occur when four or more components closely blend, which is not the case for the 
areas mapped here. A key fact that emerges from this analysis is that there is an underlying 
broad component of order 34 km/s present in each velocity regime. Verschuur & Schmelz 
(2010b) discuss the pervasive nature of this component which can clearly be seen in the 
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profiles in Figs. 4a, e & f . It is also worth noting that Gaussian analysis of all the profiles 
in the LAB survey has been carried out by Haud and described in a number of papers, see 
for example Haud & Kalberla (2007) and references therein, as well as by Nidever (private 
communication) who used an algorithm referred to in Nidever, Majewski & Burton (2008). 
Verschuur & Schmelz (2010b) have used data from both Haud (private communication) and 
Nidever (private communication) to show that results for a set of common profiles agree 
very closely and that the pervasive component does not decompose into narrow lines when 
observed with a 9 ' beam as opposed to the 36' beam of the LAB survey. Thus the Gaussian 
mapping of the HI emission from MI (and Mil, below) in the hnewidth families discussed 
here is regarded as significant. 

Figure 5 shows examples of the Gaussian decomposition for a number of directions 
toward Mil where the presence of the broad underlying component is unambiguous. After 
the Gaussian fitting was completed, a number of distinct linewidth and velocity families were 
recognized in the data. These were sorted and maps made of the total column densities for 
each prominent category. 

The average parameters for the two dominant linewidth components associated with MI 
and Mil are shown in Table 4. Column 1 gives the name of the HI feature whose profiles 
were analyzed. Column 2 gives the derived full- width, half- maximum linewidths in km/s 
averaged over the area under consideration with the 34 km/s value kept fixed. Column 3 
gives the average peak brightness temperature in K of the Gaussian components fit to the 
profiles sorted into line width family categories. Columns 4 gives the average center velocities 
in km/s with respect to the l.s.r. of these families, and Column 5 gives and average column 
densities in units of 10^^ cm~^. 

Table 4 shows that the Gaussian analysis revealed that HI emission profiles associated 
with Ml and Mil are dominated by two families of linewidths, those with this width of 
34 km/s and another of order 21 km/s wide (ranging from 18 to 26 km/s). The common 
occurrence of Gaussian components of order 34 km/s wide has been discussed by Verschuur 
(2004) and references therein and the existence of components with widths of order 21 km/s 
wide have been noted by Haud (2008) as well as the present author in the case of anomalous 
velocity HI (in preparation). It is possible that this component is due to HI at a temperature 
of 8,000 K and this will be discussed in a future paper. 



7.1. Gaussian Mapping for MI 



The dominant linewidth families found for MI have widths of 34 & 21 km/s, Table 4, 
while two other families of linewidths have average widths of 13.8 ± 2.2 and 6.3 ± 2.4 km/s. 
Figure 6a shows the ILC contours superimposed on a map of the HI column density for the 
34 km/s wide component associated with MI and Fig. 6b shows the relationship to the HI 
component with a linewidth of order 21 km/s. In this figure the southern HI and ILC peaks 
overlap nearly perfectly while the 34 km/s wide component appears to be more prominently 
associated with the northern of the two ILC peaks. Tufte, Reynolds & Haffner (1998) have 
reported excess Ha emission at five locations toward MI and these are marked in Fig. 6b 
where it is obvious, as was also noted by those authors, that the peaks in the Ha emission are 
offset from the associated bright HI peaks. When the two plots in Fig. 6 are compared with 
the map of total HI contents for MI (Fig. 4a above and more specifically Fig. 6a in Paper 
1) it is seen that it is the 21 km/s wide family of hues produces the main structure with the 
remaining contribution to the southern HI peak coming from the HI column density map 
in the two narrow components summarized in Table 4. These narrower line width values 
are consistent with the data discussed by Verschuur &: Peratt (1999). In Fig. 7 the column 
density maps for the two narrow component found in the emission profiles for MI. Together 
these associations offer a clue as to how and why the two types of emission are related, 
provided one accepts the suggestion that these associations are due to something other than 
chance. What the clues mean remains to be determined. 



7.2. A Closer Look at the HI and Soft X-ray Structure for MI 

The /LCpeak in Fig. 4a at [Ih) = (168°,67.°5) (listed as Source #9 in Table 1 of 
Paper 1) has a small X-ray structure just to its south. This is found to identically overlap 
an HI feature at positive velocities between +5 & +10 km/s at {l,h) = (169°. 67°) shown in 
Fig. 8a. In addition, an HI peak found at a velocity from —10 to —5 km/s shown in Fig. 
8b closely abuts the HI component with a 21 km/s width seen in Fig. 6b and the relative 
morphology exhibits the same axial ratio for the elongated features. These data strongly 
hint at complex interactions between HI features at different velocities interacting with one 
another to produce the HFCE observed by WMAP. Also, if the low-velocity HI feature 
associated with the slight excess of soft X-ray emission at (lb) = (169°, 67°) is significant, 
then the association with the X-ray structure for the bulk of MI is not related to the mere 
presence of high-velocity gas. 
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7.3. Maps of HI components toweirds Mil 

In the area toward Mil bounded by / = 180° & 194°, b = 62°& 68°, 225 HI profiles were 
decomposed into Gaussian components in the same way as was done for MI. Figure 9a shows 
the ILC contours overlain on a map of the HI column density in the 34 km/s wide component 
associated with Mil, and Fig. 9b shows the same contours compared to the component with 
an HI hnewidth of order 21 km/s. Tufte, Reynolds & Haffner (1998) have reported excess Ha 
emission at two locations toward MIL These are shown in Fig. 9b. Their direction #2b at 
{l,b) = (186°, 65°) lies on the peak in the HI map of the 21 km/s wide feature. The velocities 
of the HI component and Ha data are nearly identical at this position, —81.1 km/s for the 
HI and —72 or —78 km/s for each of two estimates by Tufte et al. (1998) for the Ha data. 
The underlying 34 km/s wide component at this position is centered at —82.3 km/s. 

Both of the two HI components mapped in Fig. 9 are slightly offset from the ILC peak. 
This is also true of excess soft X-ray emission associated with Mil as reported by Herbstmeier 
et al. (1995) shown in Fig. 4b. There is one bright ILC peak at {l,b) = (193.°5,63.°8) at the 
edge of the maps in Fig. 4b. It is Source 7^68 in Table 1. In Figure 10 the column density of 
a positive velocity component in the range +3 to +10 km/s is compared to the ILC contours 
and it appears that that ILC peak is associated with an excess of positive velocity HI in this 
direction. This plot also shows that positive velocity HI emission appears to encompass the 
southern half of ILC feature associated with MIL 



8. Associations in the Area Encompassing MI and Mil 

In the course of our analysis it was found that the associations between HI and the ILC 
structures are most clearly revealed in the narrow band area maps made at 2 km/s intervals. 
In Paper 1 a velocity range and interval of 10 km/s were used but in many directions this 
hides the relevant structures that cover an intrinsically smaller velocity range. 

8.1. A Striking Association south of MI 

An area of sky just south of MI includes several examples of the variety of associations 
found when small-scale features in the HI and ILC data are compared. Figure 11a shows 
the ILC contours overlain on the map of total HI content for an area bounded by I — 180° 
& 165°, and b — 45° & 65°. Many of the pairs of associated structure hsted in Paper I are 
located in this area but few are revealed in this map of total HI column density. Fig lib 
compares ILC structure and HI data integrated from —20 to —10 km/s in the same area. 
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At {l,b) — (173°, 50°) an association is visible that is listed as Source #7 in Paper 1, Table 
1 but it is barely recognized and is all but invisible in the map of total HI content in Fig. 
11a. It is shown in detail in Fig. 12a, which presents the brightness temperature for Source 
#7 in a narrow-band close-up of its area. 

A Gaussian analysis was performed on 143 profiles located every 0.°5 in latitude and 
longitude for the area bounded by / = 176° & 170°, b = 48° & 53° and again an underly- 
ing component 34 km/s wide was readily identified in most profiles for both the low- and 
intermediate-velocity HI. The most striking components that emerge after taking these into 
account is a set of narrow lines of order 3 to 5 km/s wide with a center velocity of —19 or 
—20 km/s over part of the area. The HI column density of this component was mapped 
and the result is shown in Fig. 12b with the ILC contours overlain. A near perfect overlap, 
especially as revealed in the morphological boundaries is obvious. 

If the narrow width of this component is interpreted as a kinetic temperature, it is in 
the range of 180 K to 500 K. At {l,b) = (173°, 50. °5) the narrow HI component appears to 
be a double. In Fig. 12b the column density of only one of these components was used to 
be consistent with the Gaussian solutions in its neighborhood but if the HI column density 
derived from the sum of these two components is used a map of the HI column density shown 
in Fig 12c is obtained. The bright spot at {l,b) = (172.°, 50. °5) implies that an additional HI 
structure is present there and this argues for obtaining higher angular resolution HI profiles 
for this area to better track the possible relationship to the ILC structure. In general, the 
structure seen in Fig. 12 suggests that in this direction the source of HFCE is related to 
sources associated with an enhanced region of HI emission from cold components of the gas. 

8.2. Other associations near MI 

Figure 13 illustrates the relationship between the ILC and HI structure for four other 
pairs of features in the vicinity of MI. In Fig 13a a very close association between HI at +2 
km/s and ILC source in Paper 1 is seen. In Fig. 13b the HI brightness temperature at 
+2 km/s is shown overlain with the ILC contours for Source ^^13 from Paper 1, located at 
{l,b) = (169.°5,46.°5). In that paper the positions were measured on a map made integrating 
over a 10 km/s range from —10 to km/s and the offset in position between the peaks in 
the two forms of emission was estimated at 0.°78. In the data shown in Fig. 13b, the offset 
determined from a map made by integrating over the narrower velocity band of 3 km/s is 
zero. The difference results from the presence of velocity broadening and velocity gradients 
within the structure seen in the HI. Fig. 13b also reveals another pair of associated features 
near {l,b) = (173°, 46°). This was identified as Source #24 in Paper 1 but the data shown 
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here indicate that this ILC peak is shghtly shifted with respect to the position given there. 
Figure 13c shows the HI at —136 km/s, which includes the ILC peak hsted as source # 40 
in Paper 1 at {l,h) = (167.°3,55,°7). Clearly the HI at this velocity, which peaks at {l,b) 
= (167. °3, 55. °0), is closely associated. When the difference in angular resolution between 
the two types of data are taken into account it is visually evident that these structures are 
nearly identical in shape. A Gaussian mapping of these areas was not undertaken other than 
to determine that in these directions the HI profiles showed components with fine widths of 
15 km/s for two of the peaks (from the left in Fig. 13), two overlapping components 15 and 
4 km/s wide for the third peak, and two features 34 and 20 km/s wide for the right hand 
frame. Again these hint at interacting gas masses being involved. 

Fig. 14a compares the ILC contours with the HI data at —102 km/s integrated over a 
2 km/s range. /LC Source #2 (from Paper 1) at {l,b) = (174.°4,56.°8) closely abuts an HI 
feature at {l,b) — (175. °5, 55, °0). This association between the two forms of emission for 
Source #2 is of particular interest since the boundary between the HI and the ILC peaks 
shows them to be unresolved at their interface. This is illustrated in Figure 14b, which 
displays a cross-section of the amplitudes of both the HI and ILC peaks plotted along a line 
joining the center of the two peaks. The equivalent half-widths at half-maximum intensity 
of the two features along this axis are 0.°35 and 0.°5 respectively, which compares with the 
beamwidths of 0.°6 for the HI and 1° for the ILC data. These represent a case of abutting, 
unresolved edges. The HI feature in Fig. 14a is centered at & = 55. °5 and has a total width 
of 0.°64. The ILC feature is intrinsically l.°2 wide, has a center at b =56. °5. The fine of the 
interface is thus at b = 55. °8, which is coincident with the half-maximum height of both of 
the observed features. These data suggest that the HFCE is being produced at the edge of 
the HI feature as projected on the sky. In this case the Gaussian analysis of the peak HI 
profile shows overlapping components 20 and 15 km/s wide. 



8.3. The Trio of Structures 

The basic data for a trio of features in the same general area was shown in Fig. 3. No 
attempts were made to map the Gaussian component structure for these directions, largely 
because the HI profiles are too complex for unambiguous mapping. Gaussian analysis of the 
HI peaks in Fig. 3 showed that the relevant components had widths of 7 & 4 km/s for the 
feature in (a) and 22 & 4 km/s in (b) while the IV HI emission profile in the direction of the 
HI feature in (c) was too bright to be reliably decomposed into Gaussian components. 

Overall, the data discussed in this section reinforce our contention that in order to 
carry out a comprehensive search for associations, use must be made of HI data plotted as 
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narrow channel maps spaced at intervals of at most 2 km/s and preferably with high angular 
resolution so as to untangle the profile component structure, which so far has not shown any 
simple trends. It is the HI morphology observed in small velocity intervals, and the velocity 
structure within in the profiles, that is likely to carry the most information regarding the 
nature of the complex physical interactions that appear to underlie the production of HFCE 
in interstellar space. 

9. Cross correlation studies 

Given that there appear be small angular offsets between specific ILC and HI features 
as discussed in §7 & §8, it is worth examining whether there is a statistically significant 
relationship when they are considered on a case-by-case basis. 

In §3 it was shown that there is little evidence for a direct one-to-one correspondence 
between the ILC and HI peaks over large areas of sky. Furthermore, the close associations 
discussed so far do not occur in a uniform manner because the axes connecting the members of 
paired structures arc not oriented in the same direction on the sky. Hence no comprehensive, 
large-scale statistical test for a fixed position angle over large area of sky will reveal anything 
of significance. However, a meaningful test could be to consider individual pairs of associated 
features and then shift one with respect to the other in position so as to align them even 
if they are initially slightly offset and then derive a cross-correlation coefficient to test for 
first-order significance. In this regard, small offsets between the two forms of emission would 
be the hallmark of a form of limb brightening along one face of an HI feature, where the 
neutral gas interacts with the surrounding interstellar plasma, to be discussed below. 

Table 5 shows the results of the correlation calculations for some of what are mostly 
unresolved pairs of structures discussed above and in Paper 1. Column 1 gives the source 
name and/or references used in this work and a reference to the figure number displaying 
the data. Column 2 gives the velocity of the HI at which an association is claimed while 
columns 3 & 4 give the magnitude of the cross-correlation coefficient when the two data sets 
are compared along cross sections in longitude & latitude. The values obtained using the 
data un-shifted in position are shown followed by the values obtained when the calculation 
is carried out after shifting the HI data with respect to the ILC data by the angle indicated 
in order to align the peaks. Note that if two unresolved features observed with the beam 
widths of 0.°6 and l.°0 are perfectly coincident in position the cross-correlation coefficient 
would be 0.93. Thus those cases listed in Table 5 that show a value of R of this magnitude 
are in fact highly correlated, which means that the members of each associated pair have 
virtually the same angular dimensions. It could be argued that the method of identification 
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of the closely associated pairs in the first place favors those that look closely similar. For 
example, one would tend to avoid attaching significance to finding an ILC feature of large 
angular extent next to a small HI feature, or vice versa. However, during the discovery of 
the associations it was not obvious that this was occurring. 



10. A Possible Model 

In this section it is assumed that the spatially offset associations of ILC and HI features 
is significant and will consider whether free-free emission from unresolved structures in the 
distribution of electrons in interstellar space could give rise to HFCE. Note that contrasts 
with the case of free-free emission from electron-ion pairs referred to by Hinshaw et al. 
(2007). 

10.1. Free-free Emission from Electrons 

Using the form for free-free emission as given, for example, by Nitta ct al. (1991), the 
optical depth r^j can be expressed as a function of frequency, the electron temperature, 
Te, the electron density, ng, and path length, 1, as follows: 

= 9.8 X 10"^ i/-^ T"^-^ ln{A.l x 10^° Tjv) J nl dl. (6) 
For low optical depth the observed brightness temperature as a function of u is then: 



Ts{u)^TffT,. (7) 

The expected continuum signal expected from free-free emission as defined by Eqtns. 6 
& 7 can be calculated from the angular width of the electron enhancements (or cloud) on 
the sky, 9o, and the Aspect Ratio (depth of feature relative to its width) for a given distance 
as well as the derived linear diameter or width Do- 

For a distance L in pc and expressing the scales in cm, the linear diameter is given by 

Do ^ 5.3 X 10^^ 60 L cm (8) 

If the linear depth of the electron cloud in the beam that produces the observed column 
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density is defined as D; then, 

Di = 5.3 X 10^6 OoALcm (9) 

If Nh is tfie observed HI column density in units of 10^® cm~^, then the volume density, 
uh, is given by: 

uh = 19 Nh {9o a L)'^ cm.-^ . (10) 

For a fraction, /, of the HI density in the form of electrons, the electron column density 
Ne allows the electron volume density, ng, to be derived in order to determine the emission 
measure. Here it is assumed that in general the electron "clouds" are associated with the 
adjacent HI features of known column density, although the two could be quite unrelated 
depending on how the electrons are produced. 

Using the above we find that 

rie = 18.9 / Nh {Oo A L)-^ cm-^ (11) 

Using Eqtn. (11) and substituting in Eqtns. (6) & (7), the brightness temperature as a 
function of frequency p and the observed HI column density is then given by: 

Tb{u) = 1.86 X 10^^ u-^ r-°-5 ln{4.7 x 10^° Tju) 7V| (^„ A L)-' K. (12) 

Equation 12 assumes that the source fills the beam. However, the observed brightness 
temperature will be diluted in the case of unresolved sources in the ratio of source area divided 
by the beam area. In evaluating the equation in the next section this will be considered. 

10.2. Evaluating the Equation: Another Coincidence? 

In order to evaluate Eqtn. 12 information on distance, angular scale and temperatures 
of the electron features (clouds?) is required. In the case of the most striking associations 
between HI and ILC structures shown in Paper 1 it was noted that in many cases HI at 
multiple velocities appears to be associated with a given ILC peak. Thus the HI emission 
in three velocity regimes (high, intermediate- and low-velocity) is associated. This suggests 
that all the relevant HI at high-latitudes is relatively local, within at most 200 pc of the 
Sun, if we use the canonical half-width of the galactic disk to be 100 pc and take into 



-20- 



account the latitude of the areas under consideration. This point is reinforced by the recent 
study of Verschuur & Schmelz (2010a). Furthermore, the fact that HI at distinctly different 
velocities appears to be interacting in the area where the high-frequency emission is enhanced 
suggests that the phenomenon of interaction may be involved in triggering the creation of 
electrons at an interface between the interacting HI features, or where HI features interact 
with surrounding plasma through which they travel. 

Verschuur (1991) has argued that enhancements in HI brightness seen in twisted fila- 
mentary features will be produced where the filament geometry twists into and out of the 
line-of-sight to cause the total HI column density to increase in certain directions. This de- 
fines an Aspect Ratio of A>1 for structure in filaments. However, there is no a priori reason 
to reject the possibility that A<1, which would apply in the case of sheet-like structures 
along the line of sight. 

A key parameter required to evaluate Eqtn. 12 is the electron temperature, Tg. Two 
values will be tested below. One corresponds to that of cold HI, approximately 50 K (e.g., 
Wakker et al. 1991 and a summary in Kulkarni & Heiles 1988), and 8,000 K, corresponding 
to the temperature of ionized hydrogen. 

The amplitude for free-free emission is expected to vary as z/^^ according to Eqtn. 12 
with the dependence on u encompassed by the In term small in comparison. Hence between 
the two extremes of the WMAP frequency range of 23 & 94 GHz the amphtude will decrease 
by a factor of nearly 16. However, the beam widths are 0.°88 & 0.°22 respectively and thus the 
beam dilution factor goes will increase the observed signal by a factor of 16 across the same 
frequency range, which effectively counterbalances the variation with frequency produced 
by free-free emission. Thus the resulting spectrum produced by unresolved structure in the 
beam will, to first-order, appear fiat across this frequency range. 

10.3. HFCE Produced by Unresolved Interstelleir Electron Features 

In order to evaluate the possible magnitude of this signal expected from free-free emission 
from interstellar electrons, a first-order attempt to apply observable parameters to Eqtn. 12 
taking into account beam dilution was undertaken. The HI features most closely associated 
with ILC structures for MI and Mil discussed above have a typical column density of order 
1.4 ± 0.3 X 10^° cm~^ with corresponding positive ILC average amplitude 0.16 ± 0.03 mK. 
Clearly each pair of associated structures may be at a different distance, each may have a 
different Aspect Ratio, and each may have a different temperature. However, the purpose 
of this exercise is to determine whether it is possible, for a reasonable choice of parameters. 
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to expect that HFCE at the level of 0.16 mK could be produced by free-free emission from 
clumps of interstellar electrons. 

Figure 15 plots the expected brightness temperature at both 23 and 94 GHz as a function 
of total electron column density for a number of models. Their properties were chosen to be 
illustrative only, (a) Model 1: Distance 35 pc, source diameter 6', aspect ratio 1 and — 
8,000 K. (b) Model 2: Distance 35 pc, source diameter 1', A = 0.2 and Tg = 50 K. (c) Model 
3: Distance 100 pc, source diameter 6', A = 1 and Tg =8,000 K. (d) Model 4: Distance 100 
pc, source diameter 1', A = 0.2 and Tg = 50 K. An Aspect Ratio of <1 implies a flattened 
structure along the line-of-sight. In all four plots the horizontal arrows at the center show 
the range of electron column density that would produce the observed amplitude of 0.16 mK 
and the vertical arrow shows the hmits of the brightness temperature expected at 23 & 94 
GHz at the best fit electron density. The striking fact that emerges from these examples 
is that the signals produced by free-free emission from unresolved electron structure in the 
beams at these two frequencies are closely similar and that they are relatively insensitive to 
the electron temperature. Also, the spectrum would appear nearly flat over the frequency 
range 23 to 94 GHz. The electron column densities required to produce the HFCE are of 
order 20% of the associated HI features for the MI and Mil data. 

Table 6 summarizes some of the speciflc parameter values related to these models. The 
chosen values entered into Eqtn. 12 are listed in the flrst four rows and the value for iVg 
the electron column density in the source structure required to produce the desired HFCE 
amplitude of 0.16 mK at 23 GHz is shown in the flfth row. The next four rows list the 
derived parameters for each model that are associated with the 0.16 mK signals at 23 GHz. 
The derived electron volume densities in the source for these models range from 170 to 1,400 
cm~^. These imply emission measures of order 1.6 to 12 x 10^ cm~^ pc. This may be 
compared to typical values derived from the Wisconsin H-alpha Mapper (WHAM: Haffner 
et al. 2003). Allowing for beam dilution within the WHAM beam of 1° and converting to 
Rayleighs, this would produce an observed emission measure of order 20 R, far greater than 
anything observed at high galactic latitudes, where the typical H« levels are in the range 
0.2 to 0.4 R. Thus if the WMAP structure is produced by free-free emission from unresolved 
electron density enhancements at the edges of HI features the electron temperatures cannot 
be of order 8,000 K since the corresponding Ha emission is not observed. Instead we should 
consider features where the electron temperatures may be as cold as the the HI itself, perhaps 
as low as 50K, in which case no Ha emission is expected. In that case some mechanism 
must be invoked to keep the electrons separated from the positive ions to prevent rapid 
recombination In the discussion section below we will consider what observations would be 
needed to test this hypothesis. The flnal two rows illustrate the range of electron column 
densities that are required to produce a signal of 0.16 mK at 23 and 94 GHz. The values are 
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very close, as can be seen in Fig. 15. 

11. On the nature of the possible small-scale structures 

A great deal of data relating to the small-scale structure of interstellar structure exist 
in the literature. For example, aperture synthesis observations of neutral hydrogen emission 
and absorption structure (e.g., Braun & Kanekar 2005; Diamond et al. 1989) reveal very 
small angular diameter structure to be common. Also, optical observations of absorption 
profiles produced by Nal in front of a distant stars or globular clusters (e.g., Meyer & 
Lauroesch 1999) show that the absorbing interstellar medium has structure on scales less 
than ten thousand AU across (<0.05 pc). The average angular diameter of the HI features 
found in the above study appears to be of order 1° but it is known from high resolution 
studies of high- velocity HI, for example by Schwarz & Oort (1981), Wakker (1991) and 
Wakker & Schwarz (1991), that interstellar structures can have angular dimensions as small 
as 1' in extent within larger complexes. A comprehensive discussion with many literature 
references for such "micro-structure" in the diffuse interstellar medium has been presented 
by Hartquist, Falle & Williams (2003) and in §12 we will return to relate their discussion to 
our model. 

In §10.3 it is hypothesized that unresolved concentrations of electrons in interstellar 
space could give rise to continuum emission at a level observed by WMAP summarized in 
Fig 15 and Table 6. This raises the question as to whether independent evidence exists that 
might confirm the hypothesis. The discovery by Fiedler et al. (1987) of extreme scattering 
events, or ESEs, indicate that small-scale concentration of electrons do exist. They estimated 
that electron clouds with densities of order 1000 cm~^ and a linear scale of order 0.5 AU lie 
along the lines-of-sight to several high-latitude quasars. The precise mechanism by which the 
ESEs are produced remains enigmatic, as discussed by Walker &; Wardle (1998), and from 
our point-of-view these features appear to be very much smaller than required to produce 
the necessary continuum emission observed by WMAP. 

The likely existence of the necessary small-scale structure produced as a result of MHD 
turbulence in the interstellar medium has been discussed in a series of theoretical papers, 
for example by Kowal, Lazarian & Beresnyak (2007). It is not immediately apparent how to 
relate the properties predicted by such turbulent models with the data discussed here but 
that could provide a fruitful avenue for future research 

A recent direct measurement of the electron density along a path length at a high 
galactic latitude is found in Howk, Sembach & Savage (2003) who derived a column density 
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of between 1.3 & 7.9 x lO^^^ cm-^ for a 10 kpc path toward the star vZ1128 at 6 = 78.°7. 

In part they used WHAM data. The values derived in §10.3 and summarized in Table 6 and 
Fig. 15 lie within this range. However, to reconcile these two facts implies that the bulk of 
the electron column density found by Howk et al. (2003) is confined to very small-diameter, 
high density clumps along the line-of-sight. 

To add to the circumstantial evidence in favor of deciding that wc are dealing with more 
than just odd coincidences is the fact that the properties of the structures hypothesized to 
account for the HFCE as summarized in Fig. 15 and Table 6 is that the implied volume 
densities of electrons in such features (200 to 1,400 cm~^) is in the range of the microstruc- 
turcs discussed by Hartquist, Falle & Williams (2003) said to be of order 100 AU in extent 
(equivalent to 5 x 10""^ pc. However, this is much smaller than the values pertaining to 
our model estimates, 0.03 pc at 35 pc distant and 0.01 pc at 100 pc. Since the Hartquist, 
Falle & Wihams (2003) discussion concerns neutral hydrogen structure, the question then 
becomes whether electron structures with the sizes and densities implied above also exist in 
the diffuse interstellar medium. Falle & Hartquist (2002) have argued that in a cold plasma 
a variety of slow waves can introduce large density perturbations under suitable conditions 
and in this regard we note that the diffuse interstellar medium is a cold plasma. The issue 
then becomes one of recognizing whether the HI microstructures, when they become ionized, 
will produce regional enhancements of electron density of the same order as is required to 
produce the HFCE according to Eqtn. 12. Clearly, similarities between the parameters 
implied by models such as summarized in Table 6 and the discussions of Hartquist, Falle & 
Wiliams (2003) and Falle & Hartquist (2002) deserve further consideration. 

All of these relationships can individually be described as being due to chance coinci- 
dences, which would imply that they have no physical significance. Taken together, however, 
they raise tantalizing questions. 

12. Discussion 

This paper has presented what appear to be several intriguing coincidences, (a) Small- 
scale structure observed by WMAP as summarized in the ILC map appears to be associated 
with similar structure found in the distribution of interstellar H emission when mapped in 
a small velocity band of order 3 km/s wide and separated by 2 km/s in center velocity, (b) 
The angular separation between members of associated pairs is usually of order 0.°8. (c) 
In several cases examined closely the associations appear even more dramatic when the HI 
column densities in specific line width families are mapped, (d) The angular structure of the 
two types of features are closely similar showing a cross- correlation coefficient of order unity 
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when allowance is made for beam widths and small offsets in angle, (e) In the majority of 
cases where an association has been noted, HI at more than one velocity may be involved, (f ) 
Free-free emission from interstellar electrons clumped in unresolved structures could produce 
signals of the same intensity as the high-frequency continuum emission observed by WMAP. 

Given the pervasive existence of complex HI structure in the Galactic foreground and the 
equally pervasive nature of the small-scale structure observed by WMAP, it is important to 
determine whether confusion can be created in the WMAP data by previously unrecognized 
foreground signals associated with interstellar structure. The WMAP team (e.g., Hinshaw 
et al. 2007) went to great lengths to remove contributions from well-understood sources of 
galactic radio frequency radiation, such as are observed from the galactic disk, the spurs of 
radio continuum emission, and previously mapped interstellar dust. There was no a priori 
reason for them to expect that HFCE could be produced through processes occurring in 
otherwise normal and relatively dust-free interstellar space. At the high galactic latitudes 
to which we have confined our study, very little if any cirrus dust is associated with the 
HI structures in question, although evidence for weak dust emission in HI complex M has 
recently been reported by Peek et al (2009). 

The overall pattern found in the associations found to date is that small angular offsets 
exist between paired structures. This is not unexpected, given that there is no a priori 
reason to expect that the two types of matter (neutral hydrogen and, for example, clouds 
of enhanced electron density that may give rise to the HFCE observed by WMAP), will 
coexist in identical volumes of space. More likely is the situation in which the HFCE could 
be generated at the interfaces between interacting HI features, or at the interface between 
such features moving with respect to surrounding interstellar plasma. This would produce 
a form of limb brightening observed along the region of the interface. Depending on the 
orientation of the axis describing the relative motions of the interacting gas masses with 
respect to the hne-of-sight, the associations between HI and HFCE structure may appear 
coincident in position or will more generally appear slightly offset in angle projected on the 
sky. 

If the reader is wiUing to consider that the above sections do not merely hst a string of 
odd coincidences, it may also be necessary to recognize the role of magnetic fields in con- 
centrating electrons through a pinch mechanism or to examine how electrons enhancements 
could be produced in the first place. For example, the poorly understood process discussed 
by Peratt & Verschuur (2000) may be involved in which a plasma phenomenon can trigger 
ionization when neutral material streams into a plasma. Also the role of weak dust emission 
found toward HI complex M by Peek et al. (2009) may need to be considered and such dust, 
in turn, will have a bearing on any discussion of the cooling and recombination of electrons. 
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no matter how they are produced. It is also intriguing to recognize that the volume density 
of electrons implied by the above models is of the order found in planetary nebulae but that 
the derived diameters found for these models is of order of 1 to 10% of the diameters of 
such nebulae. This would imply the existence of a category of structures smaller than the 
canonical planetary nebulae. 

Whether or not unresolved structures in interstellar electron distribution accounts for 
the high-frequency continuum emission can be tested by using very different primary beam 
widths at the same frequencies used by WMAP. Note that the Boomerang experiments (Jones 
et al. 2004) observed the small-scale structure at a frequencies of 145, 245 & 345 GHz with 
effective beam width from 11. '5 to 9.'1. If there exist structures in the electron distribution 
of order of an arcminute in extent the Boomerang data might cast additional light on this 
model. However, when the Boomerang beam widths and frequencies are entered in Eqtn. 12 
the curves in Fig. 15 are merely shifted to slightly higher electron column densities by only a 
factor of 1.3 to obtain the same average signal of 0.16 mK. Given signal-to-noise limitations, 
this may not be sufficient to test the validity of the model suggested here. 

Based on the derived values of the average electron density along the hne-of-sight, the 
models involving the source distances of order 35 to 100 pc share a common property. Using 
Eqtn. 12 it turns out that more distant features with the same internal total electron column 
density produce much weaker signals in a given beam width. For example, for Models 3 & 
4 at distances of 500 or 1,000 pc the expected brightness of the HFCE will be of order 0.03 
& 0.01 mK respectively at both 23 & 94 GHz and would be below the present threshold 
of detection. Thus even in the galactic disk foreground sources consisting of patches of 
enhanced electron density will dominate the structure observed by WMAP. 

A large-scale comparison meant to detect direct point-to-point correlations between ILC 
and galactic HI structures found in the narrow velocity bands shows no clear effect (§3), but 
then there is no obvious evidence that a point-to-point relationship is common in the data. 
It is true that foreground galactic HI features exhibit angular scales similar to those observed 
by WMAP as epitomized in the ILC map and, of course, it is well-known that the galactic HI 
emission covers all of the sky. Therefore near coincidences with HI features will be expected 
for structure in the ILC map. Thus the challenge is one again to seek independent evidence 
that could either corroborate or negate the claims made above. 

In two instances the relationship to excess soft X-ray structure reinforces the argument 
that in interstellar space there may indeed exist physical processes capable of producing 
not only some of the HFCE observed by WMAP but also excess, albeit weak, soft X-ray 
emission. In one case, MI, the result of Gaussian analysis showed that it is HI in a specific 
line width regime that is most closely associated with the ILC structure. 
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It should be stressed that the above discussion of close associations has been confined 
to structures found in an area of sky where the WMAP data required very small if any 
corrections for possible contributions from dust or synchrotron radiation as can be seen in 
the analysis of Hinshaw et al. (2007), in particular their Fig. 5. If one is prepared to 
consider that what has been found above is not due to chance associations, it will clearly be 
worth while to compare the HI data with the structures seen in the individual, un-smoothed 
data obtained by WMAP at its five frequency bands, in particular at 94 GHz, the highest 
resolution data set available. 

What has been shown above is that any attempts to confirm or negate the significance of 
the apparent associations between HI and ILC structures cannot be accomplished by simply 
comparing the ILC data with total HI column densities or even the distribution found by 
integrating over a 10 km/s velocity range as was done in Paper 1, or by Land & Slosar 
(2007). Instead, a search for closely offset associations needs to make use of HI area maps 
integrated over, say, 3 km/s and examining them at 2 km/s intervals as was done here. (Note 
that HI at 100 K will produce a line width of 2.2 km/s, which sets a practical limit to the 
resolution required.) Furthermore, such maps must be studied over the entire velocity range 
over which HI is found in any given area of sky, a range can be as large as 300 km/s, in 
order to determine at what velocity the association of a given ILC peak with an HI structure 
appears. After all, if, for example, the HFCE is produced where two HI features interact 
then the velocity at which an association is found will depend on the physical conditions of 
the HI in a given direction. These interactions will not occur at a specific velocity over the 
whole sky. 

The physics of interstellar structure, in particular of HI, is clearly complex and a full 
understanding of such structure as it pertains to creating patches of enhanced electron den- 
sities may require invoking the role of plasma and/or magneto-hydrodynamical phenomena 
as was done by Falle & Hartquist (2002) and Hartquist, Falle & WiUiams (2003). Perhaps 
this should also consider magnetic reconnection and/or plasma instabilities in field-aligned 
HI filaments in order to derive a comprehensive description of the processes that underly the 
production of HFCE in the direction of diffuse interstellar HI features. It is even possible 
that a relatively little known plasma phenomenon known as Marklund convection (Marklund 
1979), if it is operating in interstellar space, could account for separation of electrons from 
their parent ions in HI features. A full discussion of the many implications of the above 
discussions is way beyond the scope of this paper. 

If galactic foreground emission is responsible for at least some of the small-scale structure 
in the high-frequency continuum radiation observed by WMAP (and by implication COBE) 
it is then very interesting and perhaps a remarkable coincidence that the acoustic spectrum 
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said to describe the existence of sound waves in the early universe should be found to describe 
the HFCE data obtained by WMAP. 

In §11 it was concluded that the WHAM survey data appear to rule out the possibility 
that dense concentrations of electrons usually associated with a warm ionized medium could 
be present to produce the observed continuum emission observed by WMAP. Instead the 
observational restraints provided by the WHAM survey lead to the conclusion that the free- 
free emission that gives rise to a continuum signal as derived from Eqtn. 12 involves cold 
electrons at temperatures below which a significant fraction would be ionized to produce 
Ha emission. The most obvious test of the claims made in this paper then comes back 
to observing the HI and HFCE at higher resolutions to determine if apparent associations 
between HI and HFCE persist. The immediate step toward such a comparison will be to 
examine continuum data obtained by the Planck spacecraft and HI structures observed with 
the 100- meter Byrd Green Bank Telescope (GET), both of which have beam widths of order 
9'. However, even before then, those who have ready access to the highest resolution raw 
WMAP could determine if the observations at 94 GHz contain structure that can be related 
to HI data obtained with the GET. 



13. Conclusions 

Apparent associations between small-scale galactic HI features and structure in the 
WMAP ILC data have been presented above and in Paper 1. The cross-correlation between 
paired features is extremely high while a general point-to-point relationship expected for 
directly associated features is not found. It appears that the high-frequency continuum 
structure observed by WMAP may be produced by free-free emission from unresolved clumps 
of interstellar electrons. The typical offsets between the two forms of emission suggest 
that the continuum radiation originates at the interface between HI features that are either 
interacting with one another or with surrounding plasma through which they are moving. 
In general, the line-of-sight will intersect such features at some angle that would favor our 
observing an offset between the two types of radiation, although direct positional agreement 
has been noted in several cases. The possibility that these phenomena may be due to 
more than chance coincidence is reinforced by the fact that the amplitude of high-frequency 
continuum radiation is at a level expected for free-emission from electrons. 

In summary, the data considered so far suggest several alternative explanations for 
the apparent near associations between small-scale and highly correlated structures in the 
distribution of HI and the HFCE observed by WMAP. (a) The most obvious one is that 
the association are all due to chance: This includes 108 close associations reported here 
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and in Paper 1, and another 100 or so noted in a cursory examination of limited areas of 
the southern skies. However, the "chance association" conclusion is based on the unstated 
assumption that we know enough about interstellar processes that we can be sure that it is 
impossible for a previously unrecognized mechanism to produce HFCE in volumes of space 
where HI features are interacting with their surroundings, or with one another, (b) Some 
of the associations are due to chance and some may be real. This quickly poses the same 
dilemma as in (a). As soon as a few examples are shown to be clearly related to galactic 
phenomena, such as appears to be the case for MI and Mil and the examples given in Paper 
1, then the challenge is again to determine what physical mechanisms may be responsible 
for these associations so as to estimate what fraction of the ILC peaks may still be assigned 
a non-galactic origin, (c) The associations are significant. If one favors this interpretation, 
then new aspects of interstellar gas dynamics and physics may be revealed and if the above 
analysis is any indication then a full discussion will require a great deal more work than can 
be encompassed in the brief summary presented here. 

In order to understand the physical nature of the apparent associations, the properties 
of the relevant HI features must be determined by making use of higher angular and velocity 
resolution data. Thus some of the key associations discussed above should be studied using 
HI data obtained with the Green Bank Telescope with a 9' beam and compared to structure 
in the HFCE to be observed by the Planck spacecraft with a similar beam width. 

I am particularly grateful to Tom Dame for making available his program for extracting 
data from the LAB data cube for use in a Mac computer in order to produce area maps using 
a marvelous program called Transform, which was heartlessly abandoned by IDL. Without 
these two programs none of the work involved in making area maps presented here would have 
been possible on a finite budget in a finite lifetime. I wish to express my thanks to Wayne 
Landsman for providing me with the ILC data in rectangular coordinates. I am also very 
grateful for valuable discussions with and the continued encouragement from my wife, Joan 
Schmelz, who helped tide me through times when I believed that the task undertaken here 
was akin to tilting at windmills. It may still be. Ron Maddelena is thanked for suggesting 
the statistical test summarized in the §3 and Butler Burton and Mark Reid for their useful 
suggestions. Finally, an anonymous referee is thanked for extremely valuable suggestions. 
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Table 1 — Continued 
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Table 2. Top Ten ILC Peaks 



TVT 

JNo. 


7 

( 


7 




ILC lemp. 


TTT T ?■ 1 'J. 

HI Velocity 


7 

I 


7, 




TTT Tj- J 

HI amplitude 


A 1 rr 1 J- 

Angular onset 




( ) 


( ) 


(mK) 


(km/s) 


( ) 


( ) 


(K) 


( ) 








(A) 


f5) 




f7) 


f8) 


f9) 


99 


171.35 


-45.7 


0.306 


-18 


171 


-45.9 


28.6 


0.32 


99 


171.35 


-45.7 


0.306 


-30 


170.9 


-45.5 


4.2 


0.37 


100 


184.1 


-54.5 


0.259 


-50 


181.55 


-55 


0.7 


1.56 


100 


184.1 


-54.5 


0.259 


20 


185 


-56 


2.8 


1.59 


101 


160.2 


-58.25 


0.261 


-58 


160 


-58.45 


1.2 


0.23 


101 


160.2 


-58.25 


0.261 


-12 


161.5 


-58 


21.9 


0.73 


102 


176.6 


-48.95 


0.259 


-8 


177.1 


-49.45 


51.3 


0.60 


102 


176.6 


-48.95 


0.259 


12 


177.95 


-49.4 


63.5 


0.99 


103 


164.75 


-57.45 


0.252 


-6 


164.1 


-56.5 


34.6 


1.01 


103 


164.75 


-57.45 


0.252 


-66 


164 


-57 


1.3 


0.60 


104 


93.45 


-37.2 


0.251 


-8 


93 


-38 


65.2 


0.88 


104 


93.45 


-37.2 


0.251 


8 


94.5 


-36.9 


8.1 


0.89 


104 


93.45 


-37.2 


0.251 


14 


92.95 


-37.5 


4.1 


0.50 


104 


93.45 


-37.2 


0.251 


24 


93.5 


-36.9 


2.6 


0.30 


105 


94.3 


-40.8 


0.249 


18 


94.05 


-40 


2.4 


0.82 


105 


94.3 


-40.8 


0.249 


-14 


92.45 


-40.9 


22.6 


1.40 


106 


179.95 


-54.65 


0.248 


none 










107 


37.8 


42.15 


0.246 


-38 


37.05 


42 


2.2 


0.58 


107 


37.8 


42.15 


0.246 


-10 


37.9 


41.5 


12.6 


0.65 


108 


189.55 


-51.5 


0.244 


-6 


188 


-52.1 


58.6 


1.14 


108 


189.55 


-51.5 


0.244 


10 


189 


-51 


55.8 


0.61 



-34- 



Table 3. Average Separations between Paired ILC-EI Peaks 



Ref. 


/ range 


b range 


ILC peaks (mK) 


No. 


Pairs 


HI peaks 


Separation (°) 
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60°- 180° 


30°- 70° 


>0.100 
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83 


0.93 ± 0.55 
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50°- 70° 
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Table 4. Average Properties of Gaussians for Mil & MI 



Source name 


Linewidth 
(km/s) 


Peak Temperature 
(K) 


Center velocity 
(km/s) 


Column density 
(10^«cm-2) 


(1) 
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(5) 


MI 


34 
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Table 5. Correlation Coefficients 



Source name 



HI velocity 



R in longitude 



R in latitude 



Overlap 
(Paper 1) 



-87 km/s 0.81 unshifted 0.89 

0.94 offset by 0.°7 
-118 km/s 0.94 0.87 

0.94 offset by 0.°7 



Source #74 
(Fig. 14c) 

Source #79 
(Fig 14a) 

In between 
(Fig. 14b) 

Abutting peaks 
(Fig. 13) 

MI (Fig. 5b) 
22 km/s wide component 

MI 
(Fig. 7b) 



+ 1 to +3 km/s 0.71 unshifted 
0.92 offset by 0.°5 

-50 to 48 km/s 0.87 unshifted 
0.94 offset by 0.°5 



-40 to 38 km/s 



-102 km/s 



-110 km/s 



0.86 



0.95 offset by 
l.°l at an angle 

0.84 
0.89 offset 1° 



0.86 



0.87 



0.72 



0.98 resolved 
offset l°in long 



-10 to 5 km/s 0.97 resolved 0.94 

offset 0.°8 in lat offset 0.°8 in lat 



Note. — Due to beam width differences, R=0.93 for unresolved features is equivalent 
to R=1.00, see text. Structures are unresolved unless otherwise noted 
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Table 6. Model Parameters 



Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Chosen parameters: 










Electron temperature (K) 


8000 


50 


8000 


50 


Angular width (') 


6 


1 


6 


1 


Aspect Ratio 


1 


1 


0.2 


0.2 


Distance (pc) 


35 


35 


100 


100 


For Ne (10^^ cm-2) the 


34 


34 


26 


25 


derived parameters are: 










Tb at 23 GHz (mK) 


0.16 


0.16 


0.16 


0.16 


Tb at 94 GHz mK 


0.13 


0.11 


0.13 


0.10 


Volume density (cm~^) 


184 


1,100 


250 


1,400 


Average in l.o.s (cm~^) 


0.32 


0.32 


0.09 


0.08 


For Tb(z/) = 0.16 mK, 










the required A^'e values are: 










at 23 GHz (10^^ cm-^) 


34 


34 


26 


25 


at 94 GHz (10^^ cm'^) 


37 


42 


28 


31 
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Fig. 1. — Examples of the statistical tests for the relationship between the HI and ILC data 
for the northern hemisphere Target Area. The products derived from Eqtns. (4), (3) & (2) 
arc shown from left to right respectively. The product values were initially calculated in 
intervals of 0.01 K.mK and smoothed to allow this contour map representations. Contour 
intervals are 2:20@2, 25:50@5, 60:300@10 K.mK. The peak values inside the blank areas 
reach as much as 9000 K.mK. The asterisk in the center plot marks the case where the HI 
data are compared with the ILC data for the same region. See text for discussion of these 
plots. 



- 39 - 




21D 2D5 200 195 190 135 180 210 205 200 195 190 ISE 180 

Galactic longitude (") 




IDO 



150 



20C 



Fig. 2. — (a) The left-hand figure shows the total HI content integrated from —150 to +30 
km/s for an area encompassing HI feature Mil seen at {l,b)= (184°, 65°). The ILC contour 
levels overlain are from +0.03 mK in steps of 0.02 mK. In order to illustrate the challenges 
posed by comparing ILC data to the HI data, the right-hand figure (b) displays the HI total 
column density data in contour map form. In the identification of HI peaks associated with 
specific ILC peaks, the HI data in contour maps were used and the positions and amplitudes 
determined from that database. This figure does reveal several close associations in the area 
around {l,b) = (204°, 55°), which are shown in detail in Fig. 3. 
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Fig. 3. — This illustrates how the ILC contours around {l,b) = (204°, 55°) relate to the 
existence of three distinct peaks in the HI emission derived from data at three different 
velocities that are located in the area illustrated in Fig. 2. (a) The HI data in the velocity 
range —50 to —48 km/s are compared to the ILC contours. The HI structure at {l,b) = 
(201°,56.°6) is hardly visible at all in the map of total HI content in Fig. 2(a) but is directly 
associated with the ILC peak, catalogued as Source #79 in Table 1. (b) The HI data 
in the velocity range —40 to —38 km/s are compared to the ILC contours and reveals a 
clear association between HI and ILC structure at {l,b) = 202°,54.°6. (c) The HI data are 
integrated from +1 to +3 km/s and the contours in the two forms of emission are virtually 
perfectly aligned. This is source #74 in Table 1. (d) The data in the first three plots are 
here combined into one plot, which shows the very high degree of association between the 
two forms of emission for what is clearly a complex morphology at three distinct velocity 
regimes. 
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Fig. 4. — (a) The left-hand figure displays the total HI column density for MI from —140 
to —80 km/swith shaded pixels indicating the presence of excess soft X-ray emission at 1/4 
keV derived from the data of Herbstmeier et al. (1995). (b) The right-hand figure displays 
the total HI column density for Mil integrated from —120 to —55 km/s, which covers the 
great majority of HI emission at any velocity in this area of sky with the ILC contours from 
+0.02 mK in intervals of 0.02 mK overlain. The shaded pixels again indicate the excess soft 
X-ray emission derived from Fig. 7b of Herbstmeier et al. (1995). Note that the peaks found 
in all three forms of emission are slightly offset from one another in contrast to the case in 
(a) where they are closely aligned. 
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Fig. 5. — Eight frames sliowing typical profiles in the area toward Mil with the Gaussians 
fit to those profiles. In the top four frames no emission associated with Mil is found. In 
the lower four frames emission from Mil can be seen around —80 km/s. The Gaussian 
components indicated by thick lines correspond to the underlying features with a linewidth 
of 34 km/s. At low velocities these produce a near perfect fit to many of the observed 
emission profiles. In frame (g) the profile at {l,b) — (186°, 64°) shows the manner in which 
the 21 km/s wide component dominates the emission from MIL 
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Fig. 6. — The left-hand figure, (a) displays the HI column density for the 34 km/s wide 
HI component in the direction of HI feature MI compared to the ILC contours from +0.06 
mK in intervals of 0.02 mK. The double structure in the HI is uniformly offset from the 
double feature in the ILC data. In Paper 1 the total column density of the HI emission was 
compared with the ILC structure but in this plot it is evident that the physical properties 
of the HI associated with each of the two ILC peaks is distinctly different. This conclusion 
is reinforced in the right-hand plot, (b), showing the HI column density for components in 
the 21 km/s regime compared to the ILC contours. This component, indicative of HI at 
a temperature of about 8,000 K, is very closely associated with the main ILC peak. Also 
shown as filled circles are direction in which Tufte et al. (1998) detected Ha emission. Note 
that the peaks in all three categories are offset from one another, see text. 
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Fig. 7. — Two families of narrow linewidth components are also found toward MI and they 
are here plotted with respect to the ILC contours which are the same as for Fig. 6. (a) 
Column density map for components with linewidths between 9 & 15 km/s. (b) Column 
density map for narrow components with linewidth <9 km/s. Both these families appear to 
be related to the presence of the brighter of the two ILC peaks associated with MI. 
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Fig. 8. — Two {l,b) maps for the MI area again shown with respect to the same ILC contours 
used in the previous two figures, (a) In the left-hand figure a bright HI feature found by 
integrating between +5 and +10 km/s appears to be precisely aligned with the small patch 
of excess soft X-ray emission seen in Fig. 4a. See text, (b) The right-hand figure shows an 
HI feature found by integrating between —10 and —5 km/s with virtually the same angular 
extent as brightest ILC peak in this area is located just to the north of that peak by about 
0.°5. This pattern of finding HI features at two distinct velocities associated with a given 
ILC peak was found to be common, and this example is particularly striking. See text. 
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Fig. 9. — (a) At the left, the morphology of the column density of the 34 km/s wide HI com- 
ponent associated with Mil compared to the ILC contours, (b) At the right, the morphology 
of the column density of a family of linewidths in the 20 km/s range compared to the ILC 
contours from +0.02 mK in intervals of 0.02 mK overlain. Also shown as filled circles are 
direction in which Tufte et al. (1998) detected Ha emission. Unlike the case for the data in 
Fig. 6, here one of the Ha detection is closely aligned to the HI. 
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Fig. 10. — Gaussian analysis of the HI emission profiles in the area of Mil revealed unex- 
pected features at positive velocities between +5 to +20 km/s. Here the HI column density 
of these components is plotted with the same ILC contours as in Fig. 9 overlain. It is 
striking that the bright peak at {l,b) = (194°,63.°4) appears associated with a distinct ILC 
peak. Also a low level of emission at these velocities borders the brightest ILC peak in the 
area. 
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Fig. 11. — (a) The total HI column density for an area just below MI compared to the ILC 
contours level used before. Due to the dynamic range of this plot, few associations stand out. 
(b) The same area showing the integrated HI brightness from the -20 to -10 km/s at which 
resolution the associations begin to manifest. In Fig. 12 a close-up view of the structure 
around {l,b) = (173°,50.°5) shows a very high degree of association. 
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Fig. 12. — Three detailed views of the associated structures around {l,b) = (173°,50.°5) seen 
in Fig. 11. (a) A detailed view of the association between the HI brightness at —19 km/s 
(integrated over a 2 km/s bandwidth) compared to the same ILC contours as in previous 
figures. This is listed as source #7 in Table 1 of Paper 1. The outer boundary of the HI 
feature corresponds virtually perfectly to the /LC contours on three sides, (b) The HI column 
density found in the Gaussian analysis of a narrow linewidth component corresponding to 
the emission peak at —19 km/s is displayed with the ILC contours overlain. The high degree 
of correlation between the two forms of emission and the fact that they define the same area 
noted in (a) is even more striking in this plot, (c) Here both narrow linewidth components 
found in the direction marked by the bright peak are added together. In (b) only one of those 
components was included but in (c) the combined value of column densities for this double 
Gaussian (which cannot be unambiguously separated in the available HI data) overwhelms 
the morphology. 
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Fig. 13. — Four more examples in the area just south of MI that show clear associations 
between the ILC structure (same contour levels as in the previous figures) and HI emission 
brightness. Here the two forms of emission are either well aligned or slightly offset by up to 
0.°5. The two figure at the left ((a) and (b)) are both channel maps made at +2 km/s and 
the figure at the right (c) shows data at -136 km/s. 
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Fig. 14. — An example of closely abutting pair of features that appear to be unresolved at 

their interface, (a) The lower set of contours represent the HI brightness at —102 km/s (2 
km/s bandwidth) around {l,b) = (175°, 56°). The upper set of closed contours represent the 
ILC data, (b) This is a schematic cross-section along an axis joining the relevant HI and 
ILC peaks plotted in (a). The data were read off at small intervals along the axis marked in 
(a) and the scale may be judged from the fact that the HI and ILC peaks are located l.°2 
apart. At the interface between the two features, the edges are unresolved (see text) and 
notice the close agreement in overall morphology of the two forms of radiation. 
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Fig. 15. — Plots of the expected amplitude of free- free emission that would be observed at 
the two extremes of the WMAP frequency band from electron concentrations in interstellar 
space as a function of the column density of the electrons in units of 10^^ cm~^. (a) Model 

1: Distance 35 pc for electron temperature 8,000 K in a feature 6 ' in diameter and Aspect 
Ratio = 1. The horizontal line indicates the average value of the observed (positive) signal 
in the ILC data, (b) Model 3: Distance 100 pc for an electron temperature of 8,000 K for 
a feature 6' in diameter and Aspect Ratio = 0.2. (c) Model 2: Distance 35 pc, electron 
temperature 50K in a feature V in diameter and Aspect Ratio = 1 . (d) Model 4: Distance 
100 pc, electron temperature 50K in a feature V in diameter and Aspect Ratio = 0.2. See 
text. 



